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FOREWORD 


The  proper  statistical  analysis  of  rate  data  necess  Lates  tedious  and 
often  complex  calculations  involving  numerous  data  points.  These  manipu¬ 
lations  can  now  be  carried  out  easily  and  rapidly  oy  digital  computers. 

The  programs  presented  here  make  it  practical  to  carry  out  the  sta¬ 
tistical  analysis  of  kinetic  data  as  routine  part  of  a  laboratory  program, 
permitting  a  more  complete  evaluation  of  essential  chemical  properties.  The 
programming  was  prepared  and  utilized  in  the  investigation  of  compounds  for 
use  in  the  eye  protection  program.  The  work  was  performed  in  the  Flame  and 
Thermal  Protection  Section  of  the  Chemical  Modification  of  Textiles  Branch 
of  the  Clothing  and  Personal  Life  Support  Equipment  Laboratory  under  Evalua¬ 
tion  of  Chemicals  for  Optical  Shutters  -  Flashblindness  Protection  Project 
IJ662708DJ40-30.  Most  of  the  work  was  carried  out  in  the  period  April  to 
July  1971. 
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1 .  INTRODUCTION 


The  rates  of  chemical  reactions  and  the  factors  upon  which  they  de¬ 
pend  are  studied  quantitatively  using  the  techniques  of  chemical  kinetics, 
laboratory  studies  of  the  effects  of  concentration,  temperature,  and  pres¬ 
sure  on  reaction  rates  may  lead  to  conclusions  having  practical  applica¬ 
tions  or,  once  relationships  between  velocities  of  reactions  and  various 
controlling  ractors  have  been  determined,  more  fundamental  conclusions 
about  reaction  mechanisms  may  bo  arrived  at  by  interpretation  of  the  em¬ 
pirical  laws  observed  in  the  laooratory. 

In  our  investigation  of  the  triary lmethane  leuconitrile^  for  use  in 
nuclear  flashblindness  protection  devices,  it  is  necessary  to  determine  the 
kinetics  of  the  fading  reaction  of  the  leuconitriles.  Once  we  know  these 
rates  and  the  factors  affecting  them  a  system  which  returns  to  the  desired 
degree  of  transparency  in  an  optimum  period  of  time  can  be  designed. 

This  report  concerns  the  application  of  computer  techniques  to  prob¬ 
lems  of  chemical  kinetics,  that  is,  to  elucidate  the  mechanisms  of  simple 
irreversible  chemical  reactions.  Programs  which  are  presently  available  for 
evaluation  of  kinetic  data  are  generally  of  two  types.  The  first  type  in¬ 
cludes  those  that,  use  the  simplest  least  squares  analysis  and/or  are  ap- 
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plicable  only  to  simple  restricted  cases.  Although  some  of  these  pro¬ 

grams  use  advanced  statistical  methods,  each  of  them  is  applicable  only 
to  chemical  mechanisms  of  one  specific  type  (i.e.,  order). ^  At  the  other 
extreme  the  second  type  uses  very  sophisticated  statistical  methods  which 
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can  treat  even  data  front  complex  reactions  (i.e.,  reactions  whose  mechanisms 
involve  more  than  one  pathway  from  reactant  to  product,  and  in  some  cases 
involve  steady  state  intermediates  and  equilibria).  The  latter  type 
evaluates  the  non-linear  equations  involved  in  complex  reactions  by  iterative 
methods,  which  necessitate  the  use  of  large,  high  speed  computers. 

The  programs  presented  here  fill  the  gap  that  we  feel  exists  between 
the  two  above-mentioned  types  of  programs.  They  provide  programming  whose 
logic  is  readily  understandable  to  the  chemist  who  will  be  using  them,  and 
whose  simplicity  results  in.  reduced  machine  code,  core,  and  time.  They 
establish  the  rate  constants  and  order  for  irreversible  reactions  of  known 
stoichiometry,  such  as  the  triarylmethane  ion  reaction,  with  little  or  no 
restriction  on  reaction  order. 

To  present  the  programs  a  brief  description  of  chemical  kinetics  and 

Q 

its  terminology  is  given  as  background.  Then  the  programs  written  specifi¬ 
cally  for  our  work  with  the  fading  rates  of  UV  activated  photochromic  triaryl¬ 
methane  leuconitriles,  but  applicable  to  most  kinetic  problems,  can  be  dis¬ 
cussed. 

2.  DISCUSSION 

a.  Definitions 

The  rate  of  a  reaction  is  expressed  in  terms  of  the  amount  of  any 
reactant  or  product  of  the  reaction.  It  is  usually  expressed  as  change 
in  amount,  i.e.,  concentration  or  pressure,  per  unit  time.  Rates  in  solu¬ 
tion  are  expressed  as  moles  liter”*-  sec’*'  or  moles  cc~*  sec"’*'. 

The  rate  constant(k),  which  is  also  known  as  the  specific  rate  and 
as  the  rate  coefficient,  is  equivalent  to  reaction  rate  when  the  reactants 
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•re  present  at  unit  concentrations.  The  rate  constant  is  a  function  of  tem¬ 
perature,  pressure,  and  the  composition  of  the  reactant  mixture,  and,  there¬ 
fore,  indirectly  a  function  of  time.  In  stating  a  rate  or  rate  constant 
it  is  necessary  to  express  not  only  the  concentration  units,  but  also, 
the  substance,  the  reactant  or  product,  to  which  the  rate  refers. 

The  order  of  a  reaction  indicates  the  manner  in  which  the  reaction 
rate  varies  with  the  concentrations  of  the  reacting  substances.  The  order 
is  found  experimentally  by  determining  the  proportionality  between  the  rate 
and  the  power  of  the  concentration  of  each  reactant.  The  overall  reaction 
order  is  the  sum  of  the  powers  of  all  the  reactant  concentrations  which  in¬ 
fluence  the  rate. 

Thus,  for  a  reaction  of  the  following  stoichiometry 
mA  +  nb - »sC 

we  see  the  rate  can  be  expressed  by  the  changes  in  concentration: 

v  ~  “I  d(concA)  =  -l  d(concg)  =  1_  d(concG)  .  ■. 

m  dt  n  dt  s  dt  '  ' 

or 

v  =  k  (concA)8  (conCg )b  (concG)c  =  kA  (concA)a  kg  (r.oncB)b  kc  (concc)c  (2) 

where  k  =  rate  constant 
v  =  rate 

order  of  reaction  =  a  +  b  +  c 

This  equation  (eq.  2)  is  the  rate  law  for  this  particular  reaction.  In 
general,  the  form  of  the  rate  law  is  determined  by  experiment  and  cannot  be 
predicted  from  the  stoichiometric  equation  for  the  overall  reaction.  Usually 
reactions  occur  by  complex  mechanisms  which  lead  to  kinetic  laws  completely 
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different  from  the  stoichiometry.  However,  In  the  usual  case  where 
action  occurs  In  a  single  step,  there  Is  generally  a  relation  between  order 

and  stoichiometry . 

b.  The  Method  of  Integration 

The  rate  of  a  chemical  reaction  Is  usually  not  measured  directly, 
instead,  the  change  in  concentration  of  one  of  the  reactants  or  products  is 
determined  as  a  function  of  time.  The  most  widely  employed  method  of  de¬ 
termining  the  reaction  order  and  rate,  the  method  of  integration,  can  then  he 

applied. 

The  rate  of  .  reaction  may  be  related  to  the  concentration  of  one 

reactant  by  holding  all  other  concentrations  constant.  The  constant  k  then 

includes  these  other  concentrations,  and  equation  (2)  becomes: 

v  =  k  (conc)n  ^ ^ 

where  v  =  rate  observed 
ic  =  rate  constant 

n  =  order  (with  respect  to  one  reactant; 

If  for  ease  of  integration,  concentration  is  expressed  as  a0-x 
(where  a0  is  the  initial  concentration  and  x  is  the  amount  of  a  reactant 
consumed  in  time  t),  the  integrated  form  shown  in  Table  I  of  the  rate 


equations  for  different  orders  can  be  found. 


The  reaction  order  is  determined  in  this  method  by  comparison  of 
she  experimental  results  with  the  results  calculated  fur  the  integrated  rate 
equations  for  reactions  of  different  orders.  If,  as  determined  by  trial  and 


error,  the  kinetic  data  "fits"  one  of  these  integrated  equations,  (i.e., 
matches  a  calculated  curve  closely),  it  is  concluded  that  the  equation 
chosen  is  applicable.  The  rate  constant,  k,  is  obtained  for  this  order  re¬ 
action,  usually  by  put  zing  the  integrated  rate  equation  into  such  a  form  that 
a  linear  plot,  time  versus  some  measure  of  concentration  (see  Table  I, 
column  3),  is  obtained.  The  fact  that  this  second  plot  is  linear  is  a  fur¬ 
ther  indication  of  a  good  fit,  and  the  slope  (k)  can  be  found  by  the  method 
of  least  squares. 

The  method  of  integration  has  two  faults:  (1)  the  results  are  con¬ 
strained  to  fit  the  best  integral-order  law  of  those  considered,  even  if 
there  might  be  a  better  fit  with  some  other  order  and  (2)  the  reaction  is 
characterized  according  to  the  way  the  concentration  of  one  reactant  or  product 
varies  with  time.  This  technique  may  lead  to  misleading  conclusions  about 
the  way  the  rate  varies  with  the  concentration  of  the  other  reactants 
present  in  the  system. 

c.  The  Method  of  Differentiation 

The  differential  method  is  used  to  measure  the  true  order  and  the 
initial  rate  without  any  preconditions.  From  these  one  can  determine  the 
absolute  rate,  constant. 

The  differencial  method  is  based  on  the  fact  that  equation  (3)  can 
also  be  expressed  in  logarithmic  terms: 

log  v  =  log  k  +  n  log  c  (4) 

Fcr  a  v.urve  showing  the  variation  of  one  reactant  with  time,  the 
slope  at  any  given  time  and  concentration  is  the  reaction  velocity  at  chat 
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point.  Therefore,  a  double- logarithmic  plot  of  this  velocity  against  con¬ 
centration  should  give  a  straight  line  whose  slope  is  the  order  of  the  re¬ 
actant  and  whose  intercept  is  the  log  of  the  rate  for  that  order  reaction. 

TTie  initial  rate  (i.e.,  initial  slope)  of  a  number  of  runs  with  dif¬ 
ferent  concentrations  of  a  reactant  can  be  used  to  give  one  value  for  n, 
called  the  order  with  respect  to  concentration  (nc),  or  the  true  order.  Al¬ 
ternatively,  measurement  of  the  rate  at  different  times  (the  slopes  it  dif¬ 
ferent  concentrations)  during  one  run  can  be  used  to  give  another  value  for 
n,  called  the  order  with  respect  to  time  (nt).  These  two  orders  are  not 
always  the  same.  An  abnormally  large  difference  can  mean  the  reaction  is 
inhibited  (nc<nt)  or  autocatal/tic.  (nc>nt).  Both  of  these  effects  are 
due  to  interaction  of  products  with  the  reactants, 
d.  The  Isolation  Method 

rf  in  a  reaction  system  all  the  reactants  but  one  are  present  in 
large  excess,  the  apparent  order  of  the  reaction  will  be  the  order  with 
respect  to  the  one  reactant  which  is  "isolated."  This  isolated  reactant  is 
the  only  one  which  will  change  significantly  in  concentration  during  the 
course  of  the  reaction.  However,  as  with  all  these  techniques,  care  should 
be  taken  since  a  complex  reaction  may  undergo  mechanistic  changes  influenced 
by  concentration  changes,  and  the  results  may  be  unreliable. 
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3.  FORTRAN  IV  PROGRAMMING 


a.  The  Program  for  the  Method  of  Integration 


This  program  analyzes  the  kinetic  data  (i.e.,  the  time  vs.  a  changing 
variable)  in  four  :^eps.  The  first  step,  a  first  order  rate  determination 
by  least  square  linear  regression  analysis,  is  optional;  and  can  be  sup¬ 
pressed  by  use  of  the  proper  indicator.  This  optional  first  order-  rate 
analysis  is  available  and  is  printed  particularly  for  use  in  conjunction  with 
the  isolation  method.. 

In  the  second  step  of  the  program,  the  rates  at  various  times  for 
•11  the  orders  under  study  are  calculated  in  the  same  manner,  using  ex¬ 
pressions  such  as  those  found  in  Table  Tables  are  thereby  generated  for 

■  •  i 

the  rates  at  different  stages  of  the  reaction  assuming,  for  the  sake  of 

i  i 

trial,  that,  of  the  expressions  for  0,  1/2,  1,  3/2,  2,  and. 3  order  reactions, 
only  one  is  valid  for  the  experimen :al  data.  From  these  tables,  the  order 
whose  rates  at  various  times  give  the  lbwest  coefficient  of  variance  is 
selected  as  the  order. which  fits  the  experimental  data. 

In  the  third  step,  the  best  order  having  been  selected,  a  least 
square  analysis  is  used  to  determine  the  rate  of  the  reaction  for  the 
selected  order,  the  slope  of  the  line  being  obtained  by  plotting  time  against 
the  right  side  of  the  integrated  rate  expression  for  the  selected  order 


whe  re  r.^*  1 ) . 


(e,g.,  r-l  (a0  -  x) 
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Finally,  if  a  value  for  another  reactant  is  introduced,  the  integrated 
rate  expression,  for  the  reaction  (a  +  b  =  c)*  (a  second  order  mechanism  which 
is  first  order  in  a  and  b)  is  treated  tabularly  and  by  least  square  analysis  as 
described  above,  to  determine  this  second  order  rate. 

In  this  way  the  program  will  analyze  data  (1)  to  give  a  first  order 
rate  constant  for  one  reactant  for  use  if  desired  with  the  isolation  method; 

(2)  to  determine  which  of  the  trial  orders  gives  a  fit  and  determine  the  best 
rate  for  this  order;  (3)  to  calculate,  if  it  is  desired,  a  rate  for  a  second 
order  reaction  which  is  first  order  in  two  reactants. 


dX  =  kab  =  k(a0  -  x)  (bG  -  x),and  kt  = _ 1 

dt  bo  " 


In 


ao  (bn  -  x) 


bp  ( ao  “ 


(5) 


b.  The  Program  for  the  Method  of  Differentiation 

This  program  treats  the  kinetic  data  in  four  steps  to  determine  the 
order.  First,  a  relation  between  time  and  concentration  is  determined  either 
by  linear  regression  with  transformation  of  variables  (e.g.,  exponential, 
inverse,  or  semi logarithmic  curve  fit)  or  by  a  third  order  polynomial  re¬ 
gression.  Goodness  of  fit  is  evaluated  by  the  use  of  standard  errors,  and 
the  equation  which  gives  the  best  fit  for  the  experimental  data  is  used  for 
the  subsequent  steps. 

This  equation  is  then  used  in  the  secor*d  step  to  derive  the  slope,  or 
reaction  velocity,  at  various  times  including  the  initial  time.  Linear  least 


*Where  the  initial  concentrations  of  a  and  b  are  not  equal 
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squares  analysis  is  carried  out  on  the  logarithms  of  these  slopes  versus 
concentrations  at  various  times.  The  slope  of  this  regression  line  gives 
the  order  with  respect  to  time.  Next,  for  any  series  of  runs  which  differ 
only  in  the  concentration  of  the  reactant  of  interest,  the  initial  slopes 
and  initial  concentrations  for  the  series  are  treated  to  determine  a  value 
from  the  slope  of  the  regression  line  for  the  true  order  nc. 

Finally,  for  each  of  these  series,  nc  and  the  average  value  of  n^ 
are  compared.  If  they  differ  to  a  significant  extent  (  >  0.1),  the  re¬ 

action  is  labeled  either  autocatalytic  or  inhibited,  as  described  on  page  6 
above . 

In  this  way,  the  programs  analyze  kinetic  data  using  the  three  most 
common  methods  of  chemical  kinetics  to  determine  the  order  and  rate  of  re¬ 
action.  From  this  analysis  one  should  gain  a  better  understanding  of  the 
relation  between  rate  and  various  other  factors  of  the  reaction  under  study. 
This  knowledge  hopefully  would  lead  to  a  reaction  mechanism  which  will  be  of 
practical  value  in  connection  with  technical  problems. 

4,  PROGRAM  LIMITATIONS 

The  two  orograms  are  designed  (essentially)  to  treat  the  same  set  of 
input  data.  The  input  consists  of  a  variable  (such  as  concentration) 
measured  at  different  times.  A  maximum  of  twenty  time  and  vaiiable  data 
points  can  be  analyzed.  The  value  of  the  variable  must  be  greater  than  zero 
and  must  not  exceed  nine  significant  figures. 


INPUT,  OUTPUT,  AND  GLOSSARY  OF  TERNS 


Input  -  from  punched  cards 


I.  Number  of  Runs  Card  FORMAT (13) 


Column 

Variable 

1-3 

NRUNS* 

Run  Control  Card 

FORMAT(2I3, 6A5) 

1-3 

N* 

No.  of  Points 

4-6 

Nl* 

Indicator  1  =  ON 

0  =  OFF 

7-36 

CPD  (II),  II 

=1,6  Run  Name 

3.  Data  Cards  FOKMAT(8F20. 4) 

1-20, 21-40, 21-60, etc.  TB  .*1),  II  =  1,N 

1-20, 21-40, 41-60, etc.  AOTA(II),  II  =  1,N 

4.  Concentration  of  Second  Reactant  FORMAT  (F20.4) 

1-20  PO 

Output  -  on  printer 

For  each  set  of  data  the  output  appears  as: 

GIVEN  THE  FOLLOWING  XX  DATA  POINTS  FOR 


TIMF.  (SEC)  VARIABLE 

XXXXX  .XX  jHO  .XXXXE+XX 

*Right  justified 
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With  the  Nl  indicator  on,  the  program  output  from  the  first  order  kinetic 
subroutine  continues  with: 

ASSUMING  FIRST  ORDER  KINETICS,  THE  RATES  ARE 
XXXXX.XX  +O.XXXXE+XX 

MEAN  VALUE  +O.XXXXE+XX 

USING  STRAIGHT  LINE  CURVE  FITTING  PROGRAM  -  RATE  IS  +0.XXXXE+XX  AND  INTERCEPT 
IS  +0.  XXXXE+XX 

Tnen  at  the  top  of  the  next  page  the  main  program  continues: 

CONSIDERING  THIS  ORDER  WHICH  GIVES  THE  LOWEST  COEFFICIENT  OF  VARIANCE  AS  THE 

BEST  FIT,  IT  IS  FOUND  THAT  WITH  A  STANDARD  DEVIATION  OF  +0. XXXXE+XX  AND  A 

COEFFICIENT  OF  VARIANCE  OF  +0. XXXXE+XX  THE  DATA  FITS  A 

XX  ORDER  OF  RATE  EQUATION  GIVING  A  MEAN  RATE  OF  +0. XXXXE+XX  FROM  - 

XXXXX.XX  +0. XXXXE+XX 

APPLYING  A  LEAST  SQUARES  CURVE  FIT  TO  THIS  DATA 

THE  RATE  IS  +O.XXXXXE+XX 

If  none  of  the  integrated  rate  equations  have  a  reasonable  coefficient 
of  variance  (less  than  10)  the  preceding  section  is  replaced  by: 

THIS  DATA  IS  NOT  CONSISTENT  WITH  ANY  OF  THE  ONE  COMPOUND  RATE  EQUATIONS. 

The  output  then  continues  with  the  rates  as  calculated  from  all  the  integrated 
equations: 

TIME  ZERO  ORDER  HALF  ORDER  FIRST  ORDER  j/2  ORDER  2ND  ORDER  3RD  ORDER 

XXXXXX  +0.  XXXXE+XX  +0.  XXXXE+XX  +0. XXXXE+XX  +0.XXXXE+XX  +0. XXXXE+XX  +0.  XXXXE+XX 


If  a  value  for  B0  was  input  (BQ  rot  equal  0),  then  the  output  from  a 
second  order  (First  in  A  and  B)  gives: 

( ALCULATT.ON  OF  RATE  FOR  A  REACTION  FIRST  ORDER  IN  A  AND  B  WITH  AN  INITIAL 
B  OF  +0 -XXXXE+XX 

TINE  (SEC)  VARIABLE 

XXXXX.XX  jO.XXXXE+XX 

followed  by  the  Matbpac  FIT  SUBROUTINE  to  give  a  first  degree  polynomial 
curve  fit  for  a  value  of  the  rate  ohicb.  is  output  as: 

APPLYING  A  LEAST  SQUARES  CURVE  FIT  TO  THIS  DATA  THE  RATE  IS  +O.XXXXXE+XX 
However,  if  the  value  input  for  Ba  is  the  same  as  Aq,  the  output,  from  this 
subroutine  becomes: 

INITIAL  A  AND  B  ARE  EQUAL  —  SEE  SECOND  ORDER  RATES. 

After  the  last  run  has  been  processed,  END  OF  PROGRAM  is  typed  on  the  typewriter. 
Glossary  of  Terms 
Variable 
Ai,A2,A3,A4,A5 

APTAi’I),  I  =  1  ,N 
tine;i),  I  =  L ,N 
BO 

K(5,I),  I  =  1,N 
XX(I),  I  -  1 ,N 
CPD(I ) 

R(I),  I  -  1,N 
E 


li 


Description 

Dummy  arrays  for  MATHPAC/200 
SUBROUTINE  FIT 

Variable  data 

Time  data 

Value  of  initial  concentration;  of  B 

Value  of  rate  constants  for  orders  at  various  times 

Value  of  at  mjfE(I) 

Name  of  run. 

Second  order  rate  values  times  TIME (I) 

A  test  value  (ADTA(I )  -  EO) 

A  second  order  va* value 


ABR 


t"i  V^Ml  in  m 


Variable 


Description 


Sum  of  X  (XX)  values 

Sum  of  Values  of  TIME(I)  squared 

Sum  of  values  for  all  XX(I)  times  TIME(I) 

Sum  of  all  rates  times  TIME(I)  at  various  times 


SX  times  SX 


X  times  X  summed 


AMEAN 


Mean  Average  of  all  first  order  rates  at  the 
various  times 


(AN  x  SX2)  -  S2X 

Slope  of  linear  curve  fit 

Intercept  of  linear  curve  fit 


NRUNS 


Number  of  Runs 


Number  of  Points 

Indicator  to  carry  out  first  order  linear  curve  fit 


General  Indices 
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Variable 


Description 


First  Order  Rate  linear  Curve  Fitting  Subroutine 
Standard  Deviation  Subroutine 

Second  Order  Rate-determining  Subroutine  for  two 
Reactants 

Subroutine  to  Choose  Best  Order  by  Lowest 
Coefficient  of  Variance 


b.  Method  of  Differentiation  Program. 

Input  -  from  punched  cards 

1.  Number  of  Concentrations  Card  FORMAT  (13) 


CoLumf 


2,  Number  of  Runs  Card  FORMAT  (13) 


Variable 

NCONC* 


NRUNS* 


3,  Run  Control  Card  F0KMAT(I3,3X,6A5,9X,I5) 


NPTS* 


46-50 

4.  Data  Cards 


FORMAT(4F20,4) 


1- 20,  21- 40  ,  41 -60, etc. 
1-20,  21- 40, 41-  60, etc. 


CPD(I),  I.  =  1,6 
INDC* 


TIME(I),  I  -  1  ,NPTS 
ADTA(I),  J.  =  1  ,NPTS 


Output  -  on  printer 


As  in  the  integration  method,  the  run  and  time  versus  variable  data  are 
listed-  Then  the  four  curve  fitting  subroutines  each  print  the  attempted  curve 
fit,  the  best  fitting  equation  for  that  type,  and  except  ir.  the  case  of  the 
polynomial,  the  equation  for  the  slope.  Following  each  equation,  a  table  of  the 
input,  and  the  .alculated  re.sults  from  the  fitted  equation  are  printed  for  com¬ 
parison. 


*  Right  .justified. 


-**"r 


7?.w 


The  output  is: 

INVERSE  FUNCTION  CURVE  FIT 

THIS  FITS  CONC  =  1/(40 .XXXXE+XX  +  +0.XXXXE+XX  TIME)  AND 

SLOPE  =  -(40 .XXXXE+XX) (40 .XXXXE+XX  4  +0.XXXXE+XXTIME) 
and  the  table: 


-2 


TIME  VARIABLE 

40  .XXXXE+XX  +0.XXXE+XX. 


CALC.  VAR. 
+0. XXXXE+XX 


CALC.  SLOPE 
40  .XXXXE+XX 


SE  =  +0  .XXXXXE+XX 
EXPONENTIAL  CURVE  FIT 


(40  .XXXXF+XX) 

THE  FIT  GIVES  40 .XXXXE+XX  TIME 


AND  SLOPE  =  (40 .XXXXE+XX)  (+0 .XXXXE+XX )  TIME 
followed  by  a  table  like  above. 


SEMILOG  CURVE  FIT 

(40  .XXXXE+XX  +  +O.XXXE+XX  TIME) 
THIS  FITS  CONC  =10 


AND 


SLOPE  =  (2.303)  (CONC)  (+0 .XXXXE+XX) 
followed  by  a  table  like  above. 

POLYNOMIAL  CURVE  FIT 
THE  COEFFICIENTS  OF  THE  POLYNOMIAL  ARE 

nO .XXXXXE+XX  across  four  times  followed  by  a  table  as  above. 

THE  FIT  CLOSEST  TO  THE  DATA  IS 

If  any  of  the  £  rst  three  subroutines  give  the  lowest  standard  error 
they  are  relisted.  If  the  semi logarithmic  gives  the  best  fit,  then  the 
output  is: 

SEMTLOGARITHMIC  CUP YE  FIT 
and  continuing  in  all  cases: 
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!  i 


1 1 


1  1 
i  3 


I  1 


3  5* 

I 


■  -  — 


THE  DERIVATIVE  OF  THIS  EQUATION  GIVES 

TIMF  SLOPE  LOG  SLOPE  LOG  VARIABLE 

XXXXX,  XXXX  TO.XXXXE+XX  +0  „ XXXXE+XX  +C. XXXXE+XX 

At  this  point  a  polynomial  curve  fit  by  the  built-in  HATHPAC/200  FIT 

SUBROUTINE  outputs  coefficients,  variance  and  standard  error  for  the 

linear  regression,  fit  on  the  data.  The  program  continues  with: 

FROM  THIS  THE  ORDER  WITH  RESFECT  TO  TIME 

NT  =  +0  .  XXXXXE+XX 

WITH  K  =•  +XXXXXF.+XX 

If  the  time  order  indicator  is  on,  the  program  will  now  print  the  logarithms 
of  initial  concentration  and  slope: 

CONSIDERING  THE  XX  PREVIOUS  RUNS  TO  CALCULATE  THE  TRUE  ORDER 
VARIABLE  INITIAL  SLOPE 

HO  ,XXXXE+XX  +0  .XXXXE+XX 

Again,  the  MATHPAQ/ZOO  SUBROUTINE  does  a  first  order  fit  to  the  data 
printing  coefficients,  variance,  and  standard  error.  Then  the  output  continues 
FROM  THIS  THE  ORDER  WITH  RESPECT  TO  CONCENTRATION 
NC  =  HO  ..XXXXXFTOX 
WITH  K  —  HO  .XXXXXErXX 

If  the  value  differs  by  greater  than  0-1  from  the  average  Nt  already 
calculated  for  the  previous  runs, 

THE  DATA  INDICATES  THE  REACTION  TO  BE  AUTOC ATALYT T.C  (INHIBITED).. 

THE  TRlJF  ORDER  IS  TO..XXXXXE+XX 

After  the  last,  s .  t.  of  data  has  been  processed,  END  OF  PROGRAM  is  typed 
on  the  typewriter. 


■wwpppw—  1 1!  HH  'll  .IJ.1UJJ. .  . . . 


Variable 


Description 


NCONC  Number  of  concentrations  to  be  analyzed 

1NDNC  Indicator  for  determination  of  nc,  0  -  off,  1  -  on 

II 

IP 

General  Indices 

IQ 

L 

POLYN 

SEMILG  Curve  fitting  subroutines 

INVr 

EXPCF 

ERROR  A  subroutine  to  calculate  standard  error 
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TABLE  I 


Special  Cases  of  Integrated  Rate  Equations 


Order 

Differential  Equations 

Definite  Integral 

0 

dx 

dt 

=  k 

kt  =  x 

1/2 

dx 

dt 

=  k  (aQ  -  x)^ 

kt  =  2 

a0  (»o  ■ 

■  «>-*] 

1 

dx 

dt 

=  k  (a0  -  x) 

kt  =  In 

[«o/(ao 

-*>] 

3/2 

dx 

dt 

=  (aQ  -  x)3/2 

kt  =  2  ^ 

(a0  -  x; 

>4  -  (. 

2 

dx 

dt 

=  k  (aQ  -  x)2 

kt  =  x/ 

o 

C 0 

v-/ 

O 

l/Li 

-*0 

2 

dx 

dt 

=  k  (a0  -  x)  (b0  -  x) 

If 

II 

4J 

J»i 

In 

-  ao 

'•o  <bo 
bo  (ao 

3 

dx 

dt 

=  k  (a0  -  x)3 

kt  =  1/2 

2a0x 

-  X2 

ao  "  x>2 
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IZE 

^FORTRAN 

C  A  rate  ano  order  determining  prdoram  using  least  square  analysis 

C  OF  DATA  INPUT  CONSISTS  OF  TIM!  VS  REACTANT  CONC. 

C  THIS  PROGRAM  CAN  ALSO  CALCULATE  RATES  OF  A  REACTION  FIRST  ORDER  IN 

C  A  AND  IN  8.  THE  INITIAL  B  MUST  BE  INPUT  FOLLOWING  THE  CHANGE  IN  A 

C  DATA. 

COMMON  Al(8].A2t20|,A3t20l,A4t20l,A5t20J.INO.AVG,Rl20t 
COMMON  E 

01  MENS  I  ON  K 1 6. 20  I, XX 120  I, TIME  1 20 I,*DTAI20 I . A VG ( 6 ) ,CPD|6| 

REAL  K 

DATA  K,XX/140»0./ 

1  FORMAT (2 1 3, 6A5 1 

2  FORMAT) 4F20 ,41 

3  F0RMATUH1.19HGIVEN  THE  FOLLOWING, IX. 13, IX, 15HDATA  POINTS  F0R/4A5/ 
31 

4  FORMATl5X.9HT]MEtSEC1.9X.8HVARIA0LE/l 

5  FORMAT)4X,F8.2,10X,2il,4l 

6  FDRMATI1H0.16HTHE  DATA  FITS  A  . 1 2, IX , 45HORDER  OF  RATE  EQUATION  QIV 
1ING  A  MEAN  RATE  OF  , E 11 . 4. IX , 5HF RDM. I 

7  FORMATI/47HAPPLY1NG  A  LEAST  SQUARES  CURVE  FIT  TO  THIS  DATA//12X.12 
7MTHE  RATE  IS  , £12.51 

8  FORMAT  1 1M0«61HTHERE  ARE  TOO  FEW  DATA  POINTS  IlT  61  TO  CARRY  OUT  A 
2CURVE  FIT/! 

9  FOPMATI1NO,66hTHJS  DATA  IS  NOT  CONSISTENT  WITH  ANy  OF  THE  ONE  CPO 
6RATE  EOUATIONS/1 

12  FORMAT  1 1H0,36HLEAST  SOUARE  FIT  TO  THB  first  degree! 

14  FORMAT  I 90h  TIME  ZERO  ORDER  HALF  OROER  FIRST  ORDER  3/2 
70RDER  2N0  CHOER  THIRD  ORDSR/I 

15  F0RMATlI*,3x.eil.4,3X,Eli.4,3x.Ell.4.3X,Ell.4,3x.E11.4,3X.E:i.41 

16  FORMAT11H1I 

18  FORMAT I15H  END  OF  PROGRAM) 

READ  1.  NOUNS 
OD  10  1 0*1 , NRUNS 
READ  1.  N.Nl,ICPQ)III,tl*l,6) 

READ  2. ]T)M6l 1 1 I^I I-i.NI 
READ  2.  IAOTAIJJ),  JJ*l,N) 

READ  2. BO 

PRINT  3,N,!CPDIli).  1 1 • 1 7  6 1 
PRINT  4 
00  33  11*1, N 

33  PRINT  5,  TIME 1 1 1 ) ,  ADTA(II) 

IFIN1.EO.OI  00  TO  34 
CALL  SLCF(AOTA,TIME,N] 

PRINT  16 

34  DD  20  )*2.N 
X*ADTAU).AOTA)I) 

100  K)l,I|.AOTA|I)/TIMEm 
112  KI2,11«2..)ADTa)i)*ADTAI1!-ADTA(II**.5)/TIME|H 
110  K | 3, I|*ALOGlAOTA)l)/ADTA|IJi/TIMEltl 
132  K(4,  II*2.*)ADTA)!)**i:,51-ADTAH»*«i*.5tl/TIMEUI 
120  K)5,lj*X/(TIMEIII*ADTAUI*ADTAII)J 

130  K ( 6, 1  j«1 12.*A0TA)ll*XI«X*Xj/|2.*TIMEIl!*A0TA)iI*ADTAUJ*A0TAI|l*AD 


JTA  ( 1  )  | 

20  CONTINUE 

call  onriN.K.j) 
trilND.EO.76)  GD  to  91 
PRINT  6. 1  NO. AVG  t J I 
PRINT  5itTtMPtMMI.Ktj.MM).  MM»1»N) 
triN-7l70.8D.e0 
80  OD  90  JK-2.N 

XX)JK)=K)J,JK)*TiME)JK) 

90  CONTINUE 

CAL'-  PITtAl.N.O. 0.1.1. TIME. XX.A4.A5! 

PRINT  7,  A 1 1 2  I 
PRINT  16 
GO  TO  92 
70  PRINT  8 
GO  TO  10 

91  PRINT  9 

92  PRINT  l ^ 

PRINT  l5,|T|MEtt  I.Ktl.l  I.Ki2,  t  l.Kt3,  M.KI4,  t  ( .KJ5.I  )»K(6«II.  l»l.Nl 

trt80.ro. o.i  so  to  id 

CALL  SDRD|N,T)ME,'aOTa»BO( 
trtlADTAIll-aOl.EO.O.)  GO  To  10 
irtE.LE.O.OI  00  TO  lfl 
CALL  riTIAl.N.0.6.1.i.ttME,R.A4.A5l 
PRINT  7.  A 1 1 2 1 
10  CONTINUE 
TYPE  18 
STOP 
END 


1PDRTRAN 

subroutine  soRotN.TiME.Aon.aoi 

COMMON  Al(8],A2I20).A3I26l.A4l20).A5)20(.IND.AVa.R(20( 

COMMON  e 

DIMENSION  T IM1- 1 20 1 « AOTa |20  I 

4  rORMATt5X.9HT|MEtS2CI,9X.SHV*RI*BLE/l 

5  roRMAT|4X,r8.2,10X,Eii,4| 

225  PGRMATllH0.79HCALCULATI0N  OP  RATE  POR  *  REACTION  P.RST  ORDER  IN  A 
SAND  B  W)Tm  AN  INITIAL  8  OP  ,E10,4| 

226  PDRMATt50HtNlTlAL  A  AND  B  ARE  EQUAL  -  SEE  SECOND  ORDER  HATES//) 

227  PDRMATI1H0.17HB0  IS  LESS  THAN  X| 

PRINT  225, BD 

D*BO- ADTA 1 1 1 
IP  ID.EO.c, ISO  TO  221 
PRINT  4 
DO  220  )»2.N 
E»B0-A0TA(1)+A0TaI I  I 
tPIE.LE.O.CI  00  TO  230 

Rtt (•  (1,/DI ♦ALOGltAOTAtlJ*EI/)BO«ADTAl 1 1 ) I 
A8R»R  1 1 l/T  t ME ) 1 1 
PRINT  5,T)MEt I ) .ABR 
220  CONTINUE 
RETURN 


I 

t 

i 

: 

t 
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221  PRINT  226 
RETURN 

230  PRINT  227 
RETURN 
ENO 


SFQRTRAN 

SUBROUTINE  OFITIN.K.JJ 

COMMON  Al(8t.A2(20l1A3l20I.A4I20I,A5I20I,IND,AVG,RI20J 
COMMON  E 

DIMENSION  K(6,201.SDEVI6I»AVG(6] 

REAL  K 

DATA  SDEV/6»0 . / 

1  rORMATu02MCONSt0ERING  THE  ORotR  HHlcH  GIVES  THE  LOWEST  COEFFlclEN 
,  VARIANCE  AS  THE  BEST  PIT,  IT  IS  POUND  THAT//) 

SJtZ  STANDARD  DEVIATION  Of  ,Ell .  A, IX,  33«AND  A  COEFFICI 

RENT  OF  VARIANCE  OF  ,E11,4,20X. 131 
PRINT  1 
B«FL0AT[N-1I 
DO  30  J«l,  i 
VAR«0  i 
SUM.O, 

I ND  *  0 

DO  40  KK»2« N 
40  SUM«SUM*K[ J.KKI 

avgijj=sum/b 

00  30  KL*2 , N 

50  VAR«VAR*[AVGUI-K[J,KLI1«(AVGIJ]-Kt  J.KLI  I 
SDEVIJI«SQHTIVAR/B] 

30  CONTINUE 

VAR«ABS(S0EVI1I/aVGI1II 
00  60  LLE2>  6 

IFlABS[SOEV[LLl/AVG(LLIIiGT,VAR|  SO  VO  60 

VAR.ABSISDEVILLl/AVGILLl I 

J«LL 

60  CONTINUE 
I N0»  t  J-l I *5 
IFIIND.E0.25I  INo=30 
iriABSIVARI .GT.IO. I  IND»70 
PRINT  2,  SOEVIJIVVAR, IND 
RETURN 
ENO 


IFORTRAN 

subroutine  slcf(*dta,time,ni 
0 1  MENS  I  ON  AOTA (20 IjTlME 1 20  I 

11  FORMATI1H0.49HUSING  STRAIGHT  LINE  CURVE  FITTING  PROGRAM. RATE  it. IX 
3,E11.4,iX,17HAN0  INTERCEPT  IS  .Ell. 41 
i.  ^OBMAT(1H0.45HASSUMING  FIRST  ORDER  KINETICS,  THE  RATES  ARE  / J 

55  FORMAT(4X,F8.2,10X,E11,4( 

56  an«float(n'iiH**********/7X,10HME*N  V*LUE,5x»El1,41 
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. — - 


xs»o. 

SX2*0 1 

sx*o. 

SY»0  i 
S2X«0 . 

SYX»0 . 

PRINT  13 
00  210  I«2»N 
Y*AL0G[ADTAtl  1/ADTAU  )  l 
X»Y/TlMEl  I ) 

PRINT  55.  TIME  ( 1 1,’x 
XS»XS*X 

SX?»SX2*TIMEII  I  *T I  ME (11 
SYX»SYX*Y*TIMEII  I 
SY»SY*Y 
SXpSX*TIMEII) 

210  CONTINUE 
S2X»SX*SX 
AMEANaxS/AN 
B0T*AN*SX2*S2X 
A»(AN*SYX.SX«SYI/BOT 
B»ISX2*SY.SX*SYX)/R0T 
PRINT  56.AMEAN 
PRINT  11. A. 8 
RETURN 

end 


GIVEN  the  FOLLOWING  20  DATA  POINTS  FOR 
CVCN  1(1-4  ON  8,549X10-4  35C 


TiMFISECI  VARIABLE 

0.1010E-04 
0.9403E-05 
0.8754E-05 
0.8147E-05 
0 . 7629E-05 
0.7127E-05 
0 . 6415E-05 
0.6122E-05 
0.5699E-05 
0.5304E-05 
0.4932E-05 
0.4587E-05 
0.4275E-05 
0.3979E-05 
0,3708£-05 
0.3459E-05 
0 .3225E-05 
0.3007E-05 
0.2791E-05 
0.2600E-05 

ASSUMING  FIRST  ORDER  KINETICS,  THE  RATES  ARE 


1.00 

0  7190E-01 

2.00 

0.7171E-01 

3.00 

0 , 7176E-01 

4.00 

0 , 7024E-01 

5.00 

0.6981E-01 

4.00 

0.7060E-01 

7.00 

0.7158E-01 

«.00 

0.7158E-01 

9,00 

0 , 7161E-01 

10,00 

0.7172E-01 

11.00 

0 . 7179E-01 

12.00 

0 .716BE-01 

13.00 

0.7168E-01 

14.00 

0 .7160E-01 

15.00 

0.7146E-01 

14.00 

0.7137E-01 

17.00 

0 . 7129E-01 

18.00 

0.7147E-01 

19.00 

0.7144E-01 

MEAN  VALUE 

0 , 7138E-01 

USING  STRAIGHT  LINE  CyRvE  FITTING  PROGRAM-RATE  IS  0,7l58E'01  AND  INTERCEPT  IS  -0.1395E-02 


0.00 
1.00 
?.  no 

3.00 

4.00 

5.00 

6.00 

7.00 

8.00 

0.00 

10.00 

11,00 

12.00 

13,00 

14,00 

15.00 

14.00 

17.00 

18.00 

19.00 
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PPH 


CONSIDERING  7Hfi  ORPfS  WHICH  GIVES  THE  LOWEST  COEFFICIENT  OF  VARIANCE  AS  THE  BEST  FIT,  jT  IS  FOUND  THAT 


WITH  A  STANDARD  DEVIATION  OF  0.5416E.03  AND  A  COEFFICIENT  OF  VARIANCE  OF  0,75875*02  It 

THE  DATA  FITS  A  10  ORDER  OF  RATE  EQUATION  GIVING  A  MEAN  RATE  OF  0.7138E-01  FROM* 

0.00  O.OOOOE*OD 

1.00  0.7190E-01 

?.00  0.7171E-01 

3.00  0.71762-01 

4,00  0 . 7024E-01 

B.00  0 . 6981E-01 

6.00  0.7060E-S1 

7.00  0.7158E-01 

A. 00  0.715BE-01 

9.00  0 . 7161E-01 

10. 00  0 , 7177E-01 

11,00  0 , 7179E-01 

17.00  0.716BE-01 

13,00  0.716BE-01 

14.00  0 . 7160E-01 

15.00  0.7146E-01 

14.00  0 . 7137E-01 

17.00  0.7129E-01 

1B.00  0.7147E-01 

19.00  0.7144E-01 

APPLYING  a  LEAST  SQUARES  CURVE  FIT  TO  THIS  DATA 
THE  RATE  IS  0.71560E-01 
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cl 


IRMMWU'M'  "B'Up.J.. M«-  ■>' vyr-.--y. 


TIME 

ZERO  OROER 

HALF  OROER 

first  order 

3/2  OnDER 

2ND  OROER 

THIRD  ORDER 

0 

0  ,OOOOE*00 

o,ooooe«oo 

0  ,ooooe*oo 

0 iOOOOE  >00 

o,ooooe*oo 

QiO00OE«00 

1 

0 |9403E*05 

•0.6133E-02 

0 . 7190E-C1 

0 i 2303E+O2 

0.7378b*04 

0 i 7575E*0° 

2 

0.4377E-05 

•0,2959E*02 

0.7171E-01 

0 , 2339E+02 

0 i 7631E*04 

0 , 8135t *09 

3 

0.2716E-05 

•0 , 1903E.02 

0.7176E-01 

0 , P384E+02 

0 , 7925E+04 

0 , 8785F*09 

4 

0 i 1907E»05 

*  0 1 1381E*02 

0  1 7024E-01 

0,2373E*02 

0 i 8027E*04 

0 i 9233E*09 

5 

0  ■  1425E*05 

*  0 1 1 0  68E*  02 

0.6981E-01 

0i2399E*02 

0 ■ 8268E*04 

0 i 9B92E*C9 

6 

0 . 1102E-05 

*0 i 8573E»03 

0.7060E-01 

0 i 2474E*02 

0 . 870  OE*04 

0 ■ 1088E+10 

7 

0 i 8746E*  06 

*0|7069E*03 

U ■ 7l58E*01 

0»2559E*02 

0.9l96E*04 

0 i 1206E*10 

8 

0 (7124E*06 

-0.59680-1)3 

0 i 7l58g»01 

0 i 2607f*02 

0 .9562E+04 

0 i 1312f*10 

9 

0.5893E-06 

*  0  1 5l 18E»  03 

0.7161E-01 

0 ,2658E*02 

0.9952E*04 

0.1431E*10 

10 

0.4932E-06 

"0|4442E*03 

0 i 7172E»01 

0 i 2714E*02 

0.1038E+05 

0 i 1566E*10 

11 

0.4170E-06 

-0.3894E-03 

0 , 7179Ei01 

0 , 2769E+02 

0 , 1062E*05 

0  il7i5E*iO 

12 

0.3562E-06 

*0  i  3446E*  03 

0.7168E-01 

0,2818E*02 

0 i 1125E+05 

0 i 1872E*10 

13 

0 ,3061E*06 

•0.3069E-03 

0 .7168E-01 

0 ,2873E*02 

O.U72E*05 

0 , 2053E-10 

14 

0.2649E-06 

*0 i 2751E*03 

0.7160E-01 

0i2925E«02 

0 . 1219E*05 

0 i ?248E*10 

15 

0.2306E-06 

*0  |2480E»03 

0.7146E-01 

0,2974E*02 

0 i 1263E+05 

0,2459E*10 

16 

0.2016E-06 

*0  1 2245E*03 

0.7137E-01 

0|3028E«02 

0 .1319E+05 

0 ■ 2699E*10 

17 

0  1 1769E»06 

-O.20A0E-03 

0.7129E-01 

0i3083E*O2 

0 ■ 1374E*05 

0 i 2965E*1Q 

18 

0 i l55lE»  06 

•0.1656E.03 

0.7l47E:j)l 

0 i 3 155F*02 

0.144lE*o5 

0.3294E*10 

19 

0.1368E-06 

•0.1697E-03 

0.7144E-01 

0 i 3217E*02 

0 i 1503E*05 

0 i S635E*10 

CALCULATION  of  RATE  for  A  REACTION  first  ORDER  IN  a  and  8  WITH  AN  INITIAL  8  OF  0.8549E-07 


TfMFISEC) 

VARIABLE 

1.00 

0 ,8414E*02 

2.00 

0.8395E*0? 

3.00 

0.8404E*02 

4.00 

0 . 8229E*0? 

5.00 

0 .8181E*02 

6.00 

0 ,8276E*0? 

7.00 

0.8394E«02 

8.00 

0 . 8397E*02 

9.00 

0.840?E*0? 

10.00 

0 , 8418E*02 

11.00 

0 . 8428E*02 

12.00 

0 . 841 7E*02 

13.00 

0.842OE*02 

14.00 

0 . 8412E*0? 

15.00 

0 . 8398E*02 

16.00 

0 . 8389E*02 

17.00 

0 . 838lE*02 

18.00 

0 , 84O4E*02 

19,00 

0 . 8402E*02 

applying  a  lfast  squares  curve  fit  to  this  data 

THE  RATE  IS  0,8A170E*02 
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appendix  b 

the  program  for  the  method  of  differentiation 


o  o  oo  o 


IJC9  NOGO 
SF  'TRAN 

Common  Aii6).A2(30),A3t36l,A4(30(iA»[30l 
COMMON  LOATA(30I,LTIMEh0I.SEI5I.Xf301 

01  MENS  ION  CPO;6l,'TINEI30I,AOTAt30I,LSLOPI|30|,LAOTA|30l,SlNl30;,CO 
INCtSOJ.SUOPEISni 
REAL  LSLOPE ,  LAOTA,  LOATA#  LTIME 
EQUIVALENCE ILOATA.LAOTA I 
t  FORMAT! I3,3X.6A5,4X.F5.2. 1  3 1 

2  FORMAT  I 4F20 . 4 1 

3  FORMAT (lHlt19HGIVEN  THE  FOLLOWING, IX. 13, IX, 15H0ATA  POINTS  F0R/6A5/ 
31 

4  FORMaTI5X,9HT!MEISECI,V>X,0HVARIaBLE/I 

5  FORMAT  1 1H1,9X#33HTHE  FjT  CLOEST  TO  THE  DATA  IS  THE  I 

6  F0RMATI//38HTHE  DERIVATIVE  OF  THIS  EOUATION  BIVeS  /24X,4HTIME,UX, 
65HSL0PE#9x#9HL0U  St 0PE.6X, 12HL0G  VARIABLE/I 

7  FORMAT  1 21X  FjO,4,5X,3(Eli.4.5XI  I 

8  format i /40HFPOM  this  the  oroer  with  respect  to  time/i3x,5hnT  •  .Et 

B2,5//9X,9HWITK  K  *  ,E12,5I 

9  FORMAfllHl,15HCONSIOERiNq  THE.I4.42H  PREVIOUS  RUNS  TO  CALCULATE  TH 
»E  TRUE  0R0ER/9X,BHVARIaBLE.9X.13HINITIaL  SLOPE/] 

10  F0RMATI7X,Ell.4.BX,Eil.4I 

11  FORMAT  I /49HFR0M  THIS  THE  OROER  WITH  RESPECT  TO  C0NCENTRaTJ0N/13X,5 
1HNC  *  .E12.5//9X.9HWITH  K  •  .E12.5I 

12  FORMAT I//S2HTHE  OATA  INDICATES  THE  REACTION  TO  BE  AUTOCATaLTTIC.  I 

13  P0RMATI//46HTME  OATA  INOICATES  THE  REACTION  TO  BE  INHIBITEO.) 

14  F0RMATI//18HTHE  TRUE  OROER  IS  .E12.5I 

15  F0RMATI5X.F10.4.6X.E12.5I 

16  FORMAT(30x,25HSEMILOGARlfHMIC  CURVE  FJT//I 

17  F0RMAT!2lX,2?NTHlS  VARIABLE  IS  LI  0  ,///65H**N0TE**  IF  NC  IS  BEING 
7CALCULATED  FOR  THIS  SET,  IT  IS  INCORRECT.  I 

18  FORMATUSH  ENO  OF  PROGRAMI 

THE  DIFFERENTIAL  M£fHO0  Or  REACTION  OROIR  DETERMINATION  USING  A 
BEST  CURVE  FIT  FROM  E ITK2R  INVERSE,  EXPONENTIAL,  POLYNOMIAL, OR 
SEMILOGARITMMIC  CURVE  FjTS, 

READ  1.  NCONC 
DO  1000  I P»1 . NCONC 
REAO  l.NRUNS 
DT.O. 

00  1000  IQ-l.NRUNS  .  .  _ 

REAO  l,  NPTS.(CPO(Il,I.i,6l,TW.lNONC 
REAO  2,  I T I  ME (II,  I  »1^NPTS I 
REAO  2>  I AOT  A  I  1 1 ,  (cl.NPTSI 
PRINT  3,  NPTS.ICPOill.  1*1,6] 

IF  (NPTS.LT. 51  GO  TO  1000 
PRINT  4 

DO  33  1 1 *1 » NPTS 
IF  I  AOT A I  ( I  I .  LE  .  0 . 0  I  GO  TO  60 
33  PRINT15,  T I  ME  I  1 1 1 ,  AOTA(|lI 
CALL  INVYITIME.AOTA.NPTS.SLOPEI 
CALL  EXPCF  (TIME,  AOTA.’NPTS.  SLOPE  I 

call  polynitime.aota.npts.slopei 


30 


CALL  S£MILG(T{ME, ADTA.NPTS.SLOPEI 

PRINT  5 

L*1 

TEST'Sg (1 1 
no  30  1*2,4 

IFISEtl  !  .GT.TSST1  GO  TO  3ft 
1*1 

T£ST*SE 1 1 1 
30  CONTINUE 

tFtL.EQ.ll  CALL  !NVY1T|HE,ADTA,NPTS, SLOPE! 

IF (L,EQ,2(  CALL  EXPcFITIMEVaDTA.NPTS, SLOPE) 
tFIL.EQ.31  CALL  POLYN! tl HE, ADTA.NPTS. SLOPE  1 
1FIL.E0.4J  PRINT  16 
PRINT  6 

DO  20  t*l,NpTS 

IFtSLOPEUl.LT.O.t  SLOPE!  1 1*  t  '1. ) ‘SLOPE)  I ) 
LSLOPEt ) )«AL08i 0| SLOPE) 1 1 1 
LADTA|t]*ALOGlO]ADTAil)| 

PRINT  7,  T tHE 1 1 1 ,Sl OPE 1 1 1 .LSLOPEt I ) . LAOT A ( ( t 
20  CONTINUE 

CONC [ t  Q) *LADTA  (1 ) 

S ] N |  tQl*LSLOPE!l.) 

CALL  F]TlAl,NPTSil, 0,1, 1, LADTA.LSLOPB, A4, AS ) 
A1 1 1 1 *10 . **A1 1 1  ] 

PRINT  8,  A1 ] 2 ] , A1 ] 1 ) 

DT*DT*A1|2! 

IFnNDNC.EQ.nl  GO  TO  1000 
JFItQ.LT. 31  00  TO  1000 
PRINTS.  JO 

PRINT  10,  1CONC 1 1 1 , SiN  1 1  I ,  jil.NRUNSl 
CALL  FtTj*l,t0,l,l,_;i.'C0NC,S)N,A4,A3| 
A1|1|*10.**A1|1J 
PRINT  11, Al t 2 1 , A1 ] 1( 

D*A112l*DT/FLOATf]Ql 
lFIA8StDl.LT. 0.11  00  TO  40 
IF  ID  ,GT .  0  .  )  PRINT  12 
IFID.LT.o.  1  PRINT  13 
40  D*IA1[2I*DT/FL0AT1  lQtl/2. 

PRINT  14.  D 
DT*0 . 

GO  TO  10  00 
60  PRINT  17 

1000  CONTINUE  \\>e>^ 

TYPE  18 

STOP 

END 


SFORTRAN 

REAL  FUNCTION  ERROR  I ADTA,  X,  NPT8 1 
DIMENSION  AD7* 1 30 1 . X 1 30 1 
VAR»0 , 

DO  30  I  *1, NPTS 

so  var*var*iadtaj: jixu t)*)ADfAiit-x((() 


31 


ERROR'SQRTIVAR/FLOATiNPTSIl 

RETURN 

END 


SFORTRAN 

SUEROUT I NE  INVY IT J ME , AOTA. NPTS. SLOPE  I 
C  tNVERSE  CURVE  FIT 

COMMON  A1 1 6 ] ,A2l30l,A3i30l,A4l30],A5t30| 

COMMON  |.OATAI30l,LT!MEl3#l»Stl5].Xt30l 
DIMENSION  TlMB1301,ADTA(30l»SLOPi|30l 
REAL  L9L0PE,  laota,  ld*ta,  ltime 
EOUI VALENCE (LDATA.L AOTA I 

1  FORMAT I20HTH IS  FITS  CONC  ■  1/I.E11.4.3H  *  ,Ell.4,10H  TIME |  AND/66X 
1,2H.2/UX,10HSLOPE  ■  -t,Ell,4.2Hl  I.E11.4.3H  *  ,Ell.4»6H  TIMEI/I 

2  F0RMAT[//12X«4HTIME»12X,8H VARIABLE# 8X, 9HCALC .  V AR . 5X . UHCALC .  SLOP 
2E/I 

3  FORMATI10X.4IE11.4,5XI  I 

4  F0R“»Tt/////i6X,48HI  NVERSE  FUNCTION  CURVE  FI 
4T/I 

6  FORMAT ( //JHSE  «  ,'El2,5l 
00  10  1*1 i NPTS 
10  LDATA|Il»l./AOTAj(U 

CALL  FlTCAl,NPTS.O,!j,‘l,l»TlME,LOATA«A4,A5l 
PRINT  4 

PRINT  1.  A1 II 1 1 A1 1 2 1( A1 1 ? 1 <A1 1 1 1  * A1 1 2 | 

PRINT  2 

•  00  20  I»1,NPTS 

XII|*l./tAl(lI*AlI21*TtMBIIII 
SL0PEII)»*Alt2|/XII}**2 
20  PRINT  3.  TIMSlijiAOTAtiliXlII, SLOPE  1 1 1 
N ■NPTS*  2 

SEUJrERRORIAOTA.’X.NI  INVY 

PRINT  6,  SE 1 1 1 

RETURN 

ENO 


SFORTRAN 

SUBROUTINE  EXPCF IT  IMP, AOTA. NPTS, SLOPE | 

c  exponential  curve  fit 

COMMON  a1|6|,a2I30),a3(30],a4I30]#AS[30I 
COMMON  LOATAt 301, LTIME  l30],SE[5l  ,XI30l 
Oln-NSION  TIMEl30),AOTAJ30  I, SLOPE  1 30  J 
rEai  ld  at  a  ,  l  ADt  A 
EOl (VALENCE  (LOATA, LAOTA) 

1  format  i ////16X.39HE  xponentIal  Curve  fiti 

2  F0RMATI//12X,4HTJME,12X,  6HVAR I  ABLE,  8X,  9HCALC ,  VaR,  5X,  UHCALC,  SLOP 
2E/I 

3  FORMATI10X.4IE11, 4,5X1] 

5  FORMAT I/27X, Ell. 4/14HTHE  FIT  OIVIS  .E11.4.6H  TIME  //42X, 1MI, Ell, 4, 
55H  *1 i I/13HAND  SLOPE  ■  ( , Ell . 4. 2MJI, Ell . 4, 6H I  TIME//) 

6  FORMAT I//5HSE  *  , Ei2 <5 1 
00  50  1*1, NPTS 


UDATA ( I ) aALOQlO ( AOTA ( ]  I  I 
1 F IT  1  ME  1 1 1 . EQ . 0 . 1  TIME  111 >0,0  00  001 
LTIME 1 1  PAL0010 ITIMEI 1 1 1 
50  CONTINUE 

CALL  FITIA1,NPTS,0,0#1»1»LTIME.LDATA,AA,A5I 
PRINT  1 

A1UI«1C.*>A1(1I 

PRINT  5,  A1(?I,A1(1I,A1I2),A1(1I.A1IZI 
PRINT  2 

00  20  Pl.NPTS 

XII  l>AlIll«TIMEm**Al[2| 

SLOPE  I  II>A1UI*A1C2I*TIME( !  I  *  *  1 A 1 1 2 1  - 1 . 1 
20  PRINT  3,  TtME(II,ADTAtII,X(!l , SLOPE ( 1 1 
N>NPTS-2 

SE(2|>ERR0R(A0Ta,X,N* 

PRINT  6#  3EI2I 

IFITlMElH  .EO.O.oOoOOll  TIME  111*0.0 

RETURN 

END 


SFORTRAN 

SUBRO<  r(NE  POL“N(T!ME, AOTA, NPTS, SLOPEI 
C  POLYNOMI. l  curve  fit 

common  Ai(61,A2I30I»A3(30I»AA(30)»*5(30) 

common  LOATAI30I iLTIMEISOI |SE(51,X(30| 

01  MENS  ION  T I  ME  1 30  1 1 AOTA 130  I, SLOPE  1301 

1  FORMAT  1 38HTHE  COEFFICIENTS  OF  THE  POLYNOMIAL  ARE/6IE12 . 5.3X  I  ✓  I 

2  F0RMAT(//12X,AHTIME.12X,8HVAR1A8LE,8X,9HCALC.  VAR.5X, 11HCALC,  SLOP 
2E/I 

3  FORM AT (10X, 4 (Ell, 4,5X1 1 

4  FORMAT  I ////16X,37HP  QLYNOM|AL  CURVE  FIT/I 

5  FORMaT(1H0,30X,20HPOLYNOm1aL  CURVE  FIT/24X, 2 6H CANNOT  BE  CARRISO  OU 
3T  H1CH.I2.7h  POINTS) 

6  FORMAT  I //5HSE  *  .E12.5I 
IFINPTS.LE.4I  GO  TO  30 

CALL  F1TIA1,NPTS,0,0.3.3.TIME,AOTA,A4.A' I 
PRINT  4 

PRINT  1.  (All  1 1 .  1>1,4  I 
PRINT  2 

DO  20  I *1 , NPTS 

Y>T I  ME ( I  | 

XII  |*A1(1|*A1I2|«Y*A1I3I*Y>Y*A1I4  1>Y>>3*A1I5I*Y>>4 
SL0PEllI»lAl(2l*2.*Alf3l>Y*3.*Al[«I>Y.Y*4,»Al[5|>Y»3I 
20  PRINT  3,  T I  ME ( I  I , AOTA 1 1 ), X 1 1  I , SLOPE  1 1 1 
N«NPTS-4 

SE(3I»ERR0RIA0TA,X,N] 

PRINT  6,  SE 1 3 ) 

RETURN 

30  PRINT  5.  NPTS 
SE ( 3 1 >SE 12 1 *SE  ( 1  I 
RETURN 
END 


EXPCF 


POLYN 


33 


GIVEN  TUE  FOLLOWING  20  OA t A  POINTS  for 
CVCN  10»<  CN  8,549x10-4  35C 


IHEISECJ 

VARIABLE 

o.onji 

0.10100E-04 

1.000C 

0.94030E-05 

2,0000 

0.87540E-05 

3,0000 

O.8147DE-05 

4,0000 

0 , 76290E-C5 

5,0000 

0 ,7l270E-05 

6,0000 

0 , 66150E  *05 

7 .0000 

0.61220E-05 

8,0000 

0 , 56i990E*05 

9,0000 

0.53040E-05 

10,0000 

0 , 49320E- 05 

. 11,0000 

0.45870E-05 

12,0000 

0.42750E-05 

13.0000 

0 , 39790E *05 

14.0000 

0 , 37080E-05 

15,0000 

0.34590E-05 

16,0000  ; 

0 , 32250E-05 

17,0000 

0.30070E-05 

18.0000 

0.27910E-05 

19,0000 

0.26000E-05 

THIS  FITS 


INVERSE.  FUNCTION  CURVE  FIT 
CONC  »  1/t  0 , 7271E*05  «  0,146lE*05  TImEI  ANB 

SLOPE  «  *t  0,14686*05)!  0,727lE*05  *  0.146BE*05  T  t  HE  I 


-2 


TIwe 

VARIABLE 

CALC,  VAR 

CALC,  SLOPE 

0,000OE*00 

0 , 1010E»04 

0.1375E-04 

•0,7759E*14 

0 , 10  0  0E*0 1 

0.9403E-05 

0.1144E-04 

•0  1 1121E*15 

0 ,20  0  0E*01 

0.8754E-05 

0.9798E-05 

*0 , 15?9E*15 

0,3000t*01 

0 , 81 47E»05 

0 , 8566E-05 

*0.2nonE*l5 

0 , 40  0  0E*01 

0.7629E-05 

0.760VE-05 

*0 , 2534E*l5 

0 , 50  0  0E*01  ' 

0.7127E-05 

0,'6845E-05 

•0  >  3132E*l5 

0 , 6000E*01 

0.6615E-05 

0.6220E-05 

•0,3793E*l5 

0,70D0E*01 

0.6122E-O5 

0.5700E-05 

*0 , 451 7E*1? 

0.8000E*01 

0 ■ 5699E-05 

0 , 5260E«  05 

•0 ,53D4E*15 

0 ,9030E*01 

0,5304E*05 

0 i 4883E*05 

•0,6155E*15 

0,1000E*02 

0.4932E-0; 

0.4556E-05 

•0.70686*1? 

0,1100E*02 

0 ,4587g.o5 

0.4271E-05 

•0,8o45e*15 

0 . 120  0E*02 

0 , 4275E-05 

0 , 4  0l9E»05 

•0 . 9085E*1? 

0 ■ 130  0E*02 

0 .3979E-05 

0 ,3795E»05 

•0.101 9E*16 

0 ,1400E*02 

0.3708E-05 

0 ,3595E"05 

•0 , 113  6E*14 

0 , 150  0E»()2 

0 , 345  9E  *05 

0 ,34156-05 

•0.1?59E*16 

0,1600E*02 

0 .3225E-05 

0.3252E-05 

•0.138RE*16 

0,1700E»02 

0,300TE»05 

0.3104E-05 

•0.1523E»l6 

0,1«00E*02 

0 , 2791E-05 

0.2969E-05 

•0. 1665E*16 

0.1900E*02 

0 , 260  OE»05 

0.2845E-05 

•0,18l4E*l6 

SE  •  0.10519E-05 
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TIME 

0.1000E-05 

0.1000E*01 

0 ,2000E*01 

0,3000E*01 

0 i 4000E*01 

0 ,5000E*0l 

0.6000E*01 

0.7000E*01 

0,8000E*01 

0 ,9000E*01 

0,1000E*02 

0.110TE*02 

0.1200E*02 

0,l3oOE*02 

0,1400E*02 

0 ,15006*02 

0 . 160  0E*02 

0,1700E*02 

0.1800E*02 

0 ■ 1900E*0<? 


VtRUBLE 

0.1010E-04 

0 ,94o3E»05 

0 , 8754E  *05 

0,8147E*05 

0.7629E-05 

0.7127E-05 

0.6615E-05 

0 ,6t?2E*05 

0 .5699E-05 

0.530*E-05 

0.4932E-0? 

0 .4587E-05 

0.4275E*05 

0.3979E-05 

0 ,3708E*05 

0.3459E-05 

0 . 322SE-05 

0  •  3007E*05 

0.2791E-05 

0 .2600E-05 


CALC.  VAR 

0.9069E-09 

O.4945E-05 

0.7614E-05 

0.9801E-05 

0 ,11726*04 

0 ,1347E»04 

O.lSO’E'OA 

0.1661E-04 

0 ,l809E*04 

0,1943E*04 

0.2075E-0A 

0,220lE*04 

0 .232AE-04 

0.2443E-04 

O.255SE-04 

0  ,2670E*04 

0.2780E-04 

0.2887E-04 

0.2992E-04 

0.3094E-04 


calc,  slope 

0.5648E-03 
0.3079E-05 
0.2371E-05 
0.2035E-05 
0.1825E-05 
0.1678E-05 
0.1567E*05 
0.147SE-05 
0 ,1408E*05 
0.1344E-05 
0.1292E-05 
0.1246E-05 
0 ,1?06E*05 

0.1170E-05 
0 .11  38E"05 
0.1109E-05 
0.1082E-05 
0,1058E*05 
0 , 1035E*05 
0.1O14E-0S 


SE  *  0,1552«E*04 


POLYNOMIAL  CURVE  FIT 


THE  COEFFICIENTS 
0.10087F-04 


OF  THF  POLYNOMIAL  ARE  . 

0 , 69  B3lE -  06  0.20940E*07 


*0.?5837E*09 


time 

0,00006*00 
0,1000E*01 
0 ,2000E*01 
0,30006*01 
0,4000E*01 
0 ,5000E*01 
0 ,6000E*01 
0 ,7000E*01 
0,80006*01 
0,9000E*01 
0 ,10006*02 
0,1100E*02 
0 , 1?0  06*02 
0.13006*02 


variable 

0.1010E-04 
0.9403E-05 
0.8754E-05 
0.8147E-OS 
0 .7629E-05 
0.7127E-05 
0 . 6615E-05 
0 . 61 27E -05 
0 .5699E-05 
0 .5304E-05 
0.4932E-05 
0.4587E-05 
0.4275E-05 
0 .3979E-05 


Calc,  v aR 

0.1009E-04 
0.9410E-05 
0.8772E-05 
0.8174E-05 
0.7613E-V5 
0.7087E-05 
0.6596E-05 
0.6137E-05 
0.5709E-05 
0 ,5310E*05 
0.4940E-05 
0 ,  4596E  *0  5 
0.4277E-05 
0 , 3981E -05 


calc,  slope 

•0,6983E'06 
•0.6572e*06 
*0.61 74c*  0  6 
*0 .5794E-06 
-0.5432E-06 
•0.5083E-06 
-0.4749E-06 
-0.4431E-06 
*0.4129E*06 
*0 ,3842E»06 
-0.3970E-06 
-0.3314E-06 
-0.3074E-06 
-0.2B49E-06 
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!''*•*'*  '■->F>-^H'TI  '■  « ■-.-ffy  ^ 


0 .1400E-02 
0 i 1500E-02 

0,160nE*02 

0.1700E*02 
0  1  180  0E*02 
0  1 1900E*02 


0 .3708E-05 
0  i  3459E-05 
0 i 3225E-05 
0  .3007E-05 
0 . 2791E-05 
0.2600E-05 


0,3706E*05 
0.345JE-05 
0 .3217E-05 
0 .2998E-05 
D.2796E-05 
0 , 2607E-05 


-0.2639E-06 
•0 .2445E-06 
•0 , 2267E-06 
-0.2104E-06 
-0.1956E-06 
-0.1B24E-06 


SE  •  0.15926E-07 


SEMIUOO  CURVE  FIT 


THIS  FITS  CONC  «  10 


l *0 , 4995E*01  *  -0 . 31 0BE-0l  T  t HE  I 


AND  SLOPE  ■  12.303] (CONCI I-0.3108E-01J 


T I  PE 

VARIABLE 

CALC.  VAR 

CALC,  SLOPE 

0 ■ 0000F-00 

O.1O10E-04 

0.1011E-04 

-0.7238E-06 

0  a  1000E-01 

0.9A03E-05 

0 . 9416E-05 

•0 .6738f-06 

0,2000E*01 

0.8754E-05 

0 .87666-05 

-0.6272E-06 

0.3000E*01 

0.8147E105 

0.8161E-05 

*0 . 5839E- 06 

0 . 4  0  0  OE-01 

0 .7629E-05 

0.7597E-05 

-0.5436E-06 

0 ,  50  0  0E-01 

0.7127E-05 

0 • 7  073E-05 

-0.5061E-06 

o.6nooE*oi 

0.6615E-05 

0 . 65  B4E -  05 

-0.471JE-04 

0 .7000E-01 

0.6122E-05 

0.6130E-05 

•0 . 4386E-08 

0.8000E-01 

0.5699E-05 

0.57U6E-05 

•  0 ■ 4  083E- 06 

0  ,9000E*01 

0.5304E-05 

0.5312E-05 

•0.3B01E-08 

0.100DE-02 

O . 4932E- 05 

0.4945E-05 

-0.35S9E-06 

0 i 1100E-02 

0 .4587E-05 

0.4604E-05 

-0.3294E-08 

0 ■  120  0E*02 

0.4?7=Ei05 

0.4286E-05 

“0  •  3067E-06 

0  >  1300E*02 

0.39796-05 

0.399 0E-fl5 

-0.2855E-08 

0,140PE*02 

0.3708E-05 

0.S714E-J5 

-0.2658E-06 

0.1500E-02 

0 . 3459E-05 

0.3458E-05 

-0.2474E-06 

Jti6nofc*o2 

0 ■ 3225E-05 

0.3219E-05 

-0.23o3E-o6 

0,1700E*02 

0.3U07E-05 

0.2997E-05 

-0.2144E-06 

0.1fO0E*02 

0.279 lE-05 

0.2790E-05 

-0 • 1996E-06 

0.19006-02 

0.260  PE-05 

0.2597E-05 

-0.1858E-08 

SE  ■  0.193HE-07 
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fin- 


THF  FIT  CL06ST  TO  THE  DATA  IS  THE 


POLYNOMIAL  CURVE  fit 

THE  COEFFICIENTS  OF  THE  POLYNOMIAL  ARE 
0,1008 "c -04  *0 ,696316-06  0.20940E-07  *0,25837E*Q9 


time 

variable 

calc,  var 

CALC,  SLOPE 

0,000  Ob *00 

0.1010E-04 

0.1009E-04 

*0 ,6983E-06 

0 , 10  0  0£*01 

0.9403E.05 

0.9410E-05 

•0.6572E-06 

0,2000E*01 

0,87546*05 

0.8772E-05 

•0,61766*06 

0 ,30006*01 

0.8147E-05 

0.8174E-05 

-0.5796E-06 

0,40006*01 

0.7A29E-05 

0 . 761 SE-0  5 

•0.5432E.06 

0 ,5000E*01 

0,71276*05 

0 , 7  087E-05 

•0,50836*06 

0,6000E*01 

0 , 66156*05 

0,6596E»05 

-0.47496-06 

0,7000t*01 

0 ■ 61226*05 

0 . 6137E-05 

•0.4431E-06 

0.80016*01 

0 .5699E-05 

0.5709E-05 

•0.41296-06 

0,9000fc*01 

0.5304E-05 

0 , 5  3106*05 

•0.3842E-06 

0,10006*02 

0 , 4932E • 05 

0.4940E-05 

•0.3570E-06 

0,11006*02 

0.4587E-05 

0.4596E-05 

•0.33146*06 

0 , 12006*02 

0.4275E-05 

0.4?77E*05 

*0,30746-06 

0,13006*02 

0 . 3979E- 05 

0 , 3981E-05 

•0, 28496-06 

0,14006*02 

0,37086-05 

0 , 37  06E  *05 

•0.26396-06 

0,15006*02 

0.34596*05 

0.3452E-09 

•0.24456*06 

0 ,1600w*02 

0. 32256^05 

0.3217E-05 

•0.22676-06 

0,1700E*02 

0.3007E-09 

0.2998E-05 

•0.21046*06 

0,1800E*02 

0.2791E-05 

0 . 2796E  *0  5 

-0.1956E-O6 

0,19006*02 

0.2600E-05 

0 .2607E-05 

•0,1B24E«06 

SE  •  0.15926E-07 

the  derivative  of  this  equation  gives 


TIME 

SLOPE 

LOG  SLOPE 

LOG  VARIABLE 

o. ocon 

0, 69836*06 

*0.  >156E*0l 

•0 , 4996E*01 

1,0000 

0,65726*06 

*0,61826*01 

*0.5027E*01 

2,0000 

0,6176E*06 

•0.6209E*01 

•0,50586*01 

3.0000 

0.5796E-06 

•0. 62376*01 

•0,5089E*01 

4,0000 

0.5432E-06 

*0. 6265E*01 

•0.51186*01 

5.0000 

0,5o83E*06 

*0,6?94E*('l 

*0.5147E*01 

6,0000 

0.4749E-06 

*0  >  6323E*01 

*0.5179E*01 

7,0000 

0.4431E-06 

•0,63536*01 

*0 , 52l3E*01 

8.0000 

0 , 41296*06 

*0 . 6384E*0i 

•0.5244E*01 

9,0000 

0,38426*06 

•0 , 64156*01 

•0,5275E*01 

10,0000 

0.3570E-06 

•0,64476*01 

*0,53076*01 

11.0000 

0.3314E-06 

•0,6480E*01 

*0,5338E*01 

12,0000 

0.3074E-06 

-0.65126*01 

*0.5369e*01 

13,0000 

0,2849E*06 

*0 . 65456*01 

*0,5400E*01 

14,0000 

0,26396*0* 

•0.65796*01 

*0 ,5431E*01 

15,0000 

0 , 24456*06 

•0 . 66l2E*01 

•0,5461E*01 

16,0000 

0, 22676*06 

*0 , 6645E*01 

•0,549lE*01 

17,0000 

0,2104E*06 

•0,66776*01 

•0,55226*01 

18,0000 

0.1956E-06 

•0.67096*01 

•0.55546*01 

19,0000 

0,18246*06 

•0.6739E*01 

*0 . 5585E*oi 

i 

\ 


LEAST  SQUARES 


NO,  OBSERVATIONS  50 

POLY.  DEGREE  j 


COEPFICIEnTS 

n  -0 ,11412282E*0t 
1  0,1001 2S8lE*01 

ERR. SOD  t  0.72056S7E-03 

PROM  this  the  order  ,ITH  respect  to  time 
NT  =  0,10013E»0l 

WTTh  K  ■  0.722J9E-01 


polynomial  curve  fit 


STD. ErR  *  0.6327039E-02 


LEAST  SOM" 


ES  POLYNOMIAL 


CURVE  FIT 


NO,  OBSERVATION!, 
PDLV •  DEGREE 


COEFFICIENTS 

0  .0.1l593894E*02 

1  *0 .1098191Ve*OT 


RESIDUAL  SUMMARY 

POINT  X-VALUE 

1  -0 ,49956786E*01 

?  -0,4B099486E*OT 

3  -0 ,50217738E*01 


ERR-SQD  »  0.4157040E.02 


Y-VALUE 

*0 ,61559521E*01 
-0 ,6305700RE*01 
-0i60366965E*01 


y-calc 

-0,61076799E*01 

•  0 , 6311647lE*01 

•  0 ,  60790224E*  01 


STD-ERR  *  0.6447511E-01 


PROM  THIS  THE  ORDER  WITH  RESPECT  TO  CONCENTRATION 
NC  *  -t-0,10982E*01 

RjTH  K  *  0.25475E-11 


THE  data  INDICATES  THE  REACTION  TO  BE  INHIBITED. 


y.diff 

-0.48272178E-01 

0.59463233E-0? 

0.42325854E-01 


THE  TRUE  QRDFR  IS  -:.321B2E-01 
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